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Abstract 



C4__( The force-assisted desorption kinetics of a macromolecule from adhesive surface is studied 

theoretically, using the notion of tensile (Pincus) blobs, as well as by means of Monte-Carlo 
(MC) and Molecular Dynamics (MD) simulations. We show that the change of detached 
monomers with time is governed by a differential equation which is equivalent to the nonlinear 
porous medium equation (PME), employed widely in transport modeling of hydrogeological 

*^J | systems. Depending on the pulling force and the strength of adsorption, three kinetic regimes 

can be distinguished: (i) "trumpet" (weak adsorption and small pulling force), (ii) "stem- 
trumpet" (weak adsorption and moderate force), and (iii) "stem" (strong adsorption and large 
force). Interestingly, in all regimes the number of desorbed beads M(t), and the height of the 
first monomer (which experiences a pulling force) R(t) above the surface follow an universal 
square-root-of-time law. Consequently, the total time of detachment (t^), scales with polymer 
length TV as (r^) oc N 2 . Our main theoretical conclusions are tested and found in agreement 

iy~j ' with data from extensive MC- and MD-simulations. 

O 1 Introduction 



During the last decade the progress in single-molecule manipulation techniques such as atomic force 
spectroscopy (AFM) and optical or magnetic tweezers p] has attracted the interest of researchers to 
a new field of fascinating phenomena like DNA unzipping [2HH1II], or forced-induced detachment of 
individual polymers from adsorbing surfaces [5j. One can experimentally test elastic and adhesive 
properties as well as bond scission in polymer fibers, or in fundamental biological objects like 
proteins, nucleic acids, or molecular motors with spatial resolution in the nm range and force 
resolution in the pN range. In this way substantial progress has been achieved in the understanding 
and use of fibers, adhesives or biopolymers. Besides technological relevance, the field of single-chain 
manipulation poses also numerous questions of interest to fundamental physics. Thus, it has been 
shown, for example, that the force-induced unzipping of the double-stranded DNA is a 1-st order 
phase transition [21 El El [7] that has the same nature as the force-induced desorption transition of 
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a polymer from an adhesive substrate [SJ [HJ QUI ED] and that both phenomena can be described by 
the same kind of re-entrant phase diagram. 

While the theoretical understanding of these phase transformations in terms of equilibrium sta- 
tistical thermodynamics has meanwhile significantly improved [TT], their kinetics has been so far 
less well explored. Originally, a theoretical treatment of the problem was suggested by Sebastian [6] 
while the first computer simulations of the related phenomenon of DNA force-induced unzipping 
were reported by Marenduzzo et al. [7]. Yet a number of basic questions pertaining to the depen- 
dence of the total mean detachment time (r^) on polymer length N, pulling force /, or adhesion 
strength, e s , remain open or need validation. So far we are not aware of computer simulational 
studies which would shed light on these problems. Therefore, in order to fill the gap regarding 
force-induced desorption kinetics, we consider in this work a generic problem: the force-induced 
detachment dynamics of a single chain, placed on adhesive structureless plane. 

The theoretical method which we used is different and more general then in Ref. [6] where a 
phantom (Gaussian) chain was treated in terms of the corresponding Rouse equation. Instead, our 
consideration: (i) takes into account all excluded-volume interactions among monomers and, (ii) is 
based on the tensile-blob picture, proposed many years ago by Brochard-Wyart [12] as theoretical 
framework for the description of a steady-state driven macromolecule in solutions. This approach 
was generalized recently by Sakaue [T3l [H] for the analysis of non-equilibrium transient regimes in 
polymer translocations through a narrow pore. 

Very recently, Sebastian et al. [15] demonstrated (within the tensile blob picture) that the average 
position x(n, t) of a monomer n at time t is governed by the so-called p-Laplacian nonlinear diffusion 
equation [16]. In contrast, in the present paper we discuss force-induced desorption kinetics on the 
basis of a nonlinear diffusion equation which governs the time-dependent density distribution p(x, t) 
of monomers. We show that this equation is similar to the porous medium equation (PME) [IT] 
which has gained widespread acceptance in the theory of groundwater transport through geological 
strata [18]. The most natural way of solving the PME is based on self- similarity which implies that 
the governing equation is invariant under the proper space and time transformation [T71 [18] . Thus, 
self-similarity reduces the problem to a nonlinear ordinary differential equation which can be solved 



much more easily. 

Eventually, in order to check our predictions we performed extensive MD- and MC-simulations 
which, for the case of strong adsorption, large pulling forces and overdamped regime, show an 
excellent agreement with our theoretical predictions. Moreover, the comparison of MC and MD 
data reveals that the presence of inertial effects, typical in underdamped dynamics, significantly 
affects the observed detachment kinetics. 

The paper is organized as follows: after this short introduction, in Section [2] we present briefly 
the theoretical treatment of the problem and emphasize the different regimes pertinent to polymer 
detachment kinetics. A short description of the simulational MC- and MD models is given in 
Section [3] while in Section H] we discuss the main results provided by our computer experiments. In 
Section \5\ we suggest a plausible experiment which could validate our theoretical results. We end 
this work by a brief summary of our findings in Section 

2 Tensile-blob picture for the force-induced detachment 

2.1 Trumpet regime 

2.1.1 Beads density temporal variation is described by a nonlinear porous medium 
equation 

After equilibration on an attractive solid substrate, the adsorbed polymer chain can be pictured as 
a two dimensional string of adsorption blobs of size D oc a(e s — e c )~ u [IHL where a is the Kuhn 
segment length, and v and <fi are respectively the Flory exponent and the crossover (or, adsorption) 
exponent. The dimensionless adsorption energy e s = Eg/ksT (here and in what follows k^ denotes 
the Boltzmann constant, and T is temperature), whereas e c is the critical adsorption point (CAP). 
The external driving force / acts on the free chain end (in the x-direction) while the second chain end 
is tethered to the substrate. If the pulling force strength falls in the range 1/N V <C aj jk^T ^ 1, the 
chain tail forms a kind of "trumpet" shape, as shown schematically in Fig. [T^,. Moreover, as far as 
the adsorption blob, adjacent to "trumpet", changes its lateral position (i.e., its coordinates within 
the substrate) during desorption, the "trumpet" has not only vertical velocity (in the x-direction) 
but also a lateral one. The latter leads to a tilt of the "trumpet" axis in direction opposite to the 
lateral velocity due to Stokes friction. The tilting angle is denoted by 9 in Fig. [TJ and, as will be 



seen from our simulation results in Sec. H] almost does not change in the course of detachment. 
This substantially simplifies the theoretical consideration. 
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Figure 1: The picture of the single adsorbed chain detachment. The chain is initially equilibrated, 
giving rise to a string of adsorption blobs with a blob size D. The external driving force /, exerted 
on the first monomer, is switched on at the initial time moment. By time t the position of the 
first monomer is denoted as R(t). The shaded part of the chain displays a tensile blobs sequence 
which has a form of a "trumpet" , tilted opposite to the lateral velocity direction by angle 9. a) The 
"trumpet" regime, corresponding to weak adsorption and weak pulling force, i.e., 1/N V <C <j? < / < 
1. b) The "stem" regime, corresponding to strong adsorption and strong pulling force, 1 <C (p < f- 
c) The "Stem-trumpet" regime, representing the intermediate pulling force interval tp < 1 < /. 



Let us denote by f(x,t) the tensile force which acts in direction of the "trumpet" axis at an 
arbitrary cross-section of a tensile (Pincus) blob, placed at distance x from the substrate at time t. 
The x-component of this tensile force reads f x (x,t) = f(x,t) cos6. The blob size is then given by 
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(1) 



f(x,t) f x (x,t) 

where C = cos 9. It can be seen from Fig. [Th that the size of the tensile blob at x corresponds to 
the lateral size of the polymer at height x. The balance of forces in x-direction at this point yields 



fx{x,t) =7 / v x (x',t) 
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where v x (x, t) is the x-component of the local velocity and (p is the restoring force which is exerted 
on the moving domain by the last adsorption blob. The value of (p depends on the proximity to 
the critical point e c as well as on the chain model itself (see below). In the integral in Eq. 02]), 
dx'/£(x', t) counts the number of blobs in the interval x', x' + dx' whereas jv x (x', t) [£(V, t)/a] z ~ 2 is 
the local Stokes friction force (7 being the friction coefficient) [2D]. The dynamic exponent equals 



either z = 2 + 1/V, or z = 3, for Rouse- and Zimm dynamics respectively. Thus, taking into account 
Eq. ((!]), the equation for the blob size reads 



Ck B T 
7 P D v x (x',t) [f(x')/a] z ~ 2 <&7f (*') + ¥> ' 



&> = rx ^ ,w 18 -2, ~ • ( 3 ) 



With this equation one may represent the local velocity v(x,t) in terms of the blob size £(x,t) as 
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The local (linear) monomer density p(x,t) and the local blob size £,{x,t) are related by p(x,t) = 
[€{x,t)/a]( 1 ~ v '' u /a so that Eq. (j3j) can be written in the form 

,; f T f) = CkBTV l 9 I^1 (<\ 

x{,) 7(1-1/) C9(x,t)]"(*- 2 )/( 1 -")+ 1 ax ' lJ 

where we have introduced the notation for the dimensionless monomer density p(x,t) = ap(x,t). 
The local flux of polymer beads has a standard form J x (x, t) = v x (x, t) p(x, t), that is, 

T ( ,x uCk B T d_ 

x{ ' J ~ a 7 (l - u)[p(x,t)Y^)/(^) dx P{X) } [ ' 

With J x (x,t) from Eq. (jBJ), one may determine p(x,t) from the continuity equation dp(x,t)/dt = 
—dJ x (x : t)/dx: 

m dx v 9xy 

where we have introduced the dimensionless (tilded) independent variables x = x/a and t = t/ro 
(thereby the characteristic time is To = a 2 7(l — v) j vCk B T\ Eq. (J2J) has the form of the nonlinear 
porous medium equation (PME) [17] where the characteristic exponent m = —[v{z — 1) — 1]/(1 — v) 
has in our case a negative value. An equation similar to PME has been discussed recently by Sakaue 
et al. [2TJ in the context of polymer expansion and unfolding from a compact state. In this case 
the consideration has been based on the dynamics of concentration blobs (20] rather than tensile 
blobs like here. This kind of nonlinear partial differential (parabolic) equations is commonly used 
in the context of groundwater infiltration or seepage through porous strata [18] where the exponent 
m > 1. 



The geometric factor C = cos 9 is unessential in the subsequent consideration and we will drop 
it. In addition to Eq. (J7J), one has to fix the initial and boundary conditions (BC). Initially, there 
are no polymer segments at x > D, so that 

p(x, t = 0) = (8) 

At x = D the tensile blob is equal in size to the adsorption blob, i.e., £(x = D,t) = hsT/ip. As a 
result, in terms of the dimensionless density, one has 

p(x = D,t) = — (9) 

where D = D/a and cp = aipjksT ■ 

The BC on the moving end of the "trumpet" follows from the condition £(x = R(t),t) = ksT/f 
which corresponds to the smallest blob size, i.e., 

p(x = R(t),t) = 1 J—- (10) 

Here again the dimensionless quantities, R(t) = R(t)/a and / = af/ksT, are introduced. 

Our consideration is based on the assumption that the lateral component of the tensile force 
is very weak (because the trumpet is almost vertical), so that any force propagation along the 
adsorbed part of chain can be neglected. This is different from the case considered by Serr and 
Netz [22] where the pulling of a strongly adsorbed chain by means of an AFM has been investigated. 
In contrast to the present case where an uniform field pulls on the first bead, the AFM cantilever 
can move vertically and horizontally with respect to substrate so that the tilt angle 9 changes with 
time. Moreover, in this case the sliding of the polymer on the substrate becomes important with the 
polymer- surface friction force being much larger than the Stokes friction force for the detached part 
which is essential in our case. Therefore, the AFM can not be actually used for testing the scheme 
shown in Fig. [TJ In Sec. |5]we will discuss a plausible electrostatic tracking experiment which makes 
it possible, in principle, to test a gentle detachment even close to the adsorption critical point. 

2.1.2 Self-similar solution 

In the same manner as for the PME, the solution of Eq. (j7|) is derived by symmetry of self-similarity 
(SSI. This svmmetrv implies invariance of the solution with respect to the stretchino transformation 



: t — > r/t , x — > r] S x, p — > rf p. The requirement that Eq. ([7]) stays invariant under the stretching 
transformation connects the exponents u,m, and 5 by the following relation [17] 

cu(l-m) + 26=l (11) 

On the other hand, this requirement also fixes the self-similar form of the solution as 

p(x,t) = t u W(u), u = (x-D)/t 5 . (12) 

In Eq. (TT2I) . W(u) is a scaling function and we took into account that at t = x = D. On the 
other hand, the boundary condition given by Eq. (Q is compatible with the SS-form Eq. (TT2l) only 
if oj — 0. Then, as a consequence of Eq. ffTTl) . 5 = 1/2. Substituting Eq. (fT2"|) into ([7]) yields a 
nonlinear ordinary differential equation for the scaling function W(u), i.e., 

d( W m-ld w \ + ud w = Q 

du \ du J 2 du 

The BC given by Eq. ([9]) reads 

mu = oj = ^ (u) 

Assume now that the moving chain end follows the scaling relation R(t) = D + u^t x where x is 
an exponent and the factor v^ will be fixed below. This scaling ansatz is compatible with the BC 
given by Eq. ( flOl) only if \ = 1/2. This imposes the boundary condition 

W{u = v d ) = -J— (15) 

f{l-v)/v 

Finally, the initial condition Eq. ([8]) becomes 

W(u = oo) = (16) 

The nonlinear Eq. (1131) may be solved numerically. However, in case the pulling force / exceeds 
the restoring force only by a small amount, one can show (see below) that Vd <^C 1, i.e., the interval 
of argument u variation, where W(u) is nonzero, is very narrow. As a result, < u < v^ "C 1 and 
one can linearize Eq. ( TT3l (by making use Eq. ( Tl4j) ) to get 

d ii d 

V z - 2 ^W(u) + --W(u) = (17) 



where for small arguments W(u) m 1 ~ p( z 2 \ The solution of Eq. ( I17p . subject to both conditions 
Eq. (TT4l) and Eq. (fl6l) . is given by 
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where erf(x) = (2/y/n) J* exp(—t 2 )dt is the error function. 

The chain end velocity amplitude, v^, can be determined now from the condition Eq. ([15 
If the pulling force / is only slightly larger than the restoring force tp, one can expand W(u) '< 
[1 - uK-n 1 ' 2 ^- 2 ^ 2 )}/^ 1 -^^. Thus, one obtains 



V d ~ y/lTp 



-(z-2)/2 



(l-v)/v 



The equation of motion for the moving end then reads 



R(t) = D + ^p 



(*-2)/2 



(l-u)/u 
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(19) 
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The variation of the number of desorbed monomers with time, M(t), is given by the integral 

R(t) v d 

M(t) = / p(x,t)dx = Vt W(u)du 
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Eventually, the total time for polymer detachment, (r d ), reads 

r iV 2 



(r.) 
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Combining Eqs. (12 ip and ( 1221) . one obtains the following dynamical scaling law 

M(t) = Ni 



t 



(r d ) 



(23) 



Eventually, we would like to stress that in the foregoing consideration one neglects the thermal 
fluctuations so that the necessary condition for desorption requires that / > (p. On the other hand, 
the "trumpet" formation is only possible when 1/N U <C / < 1. Therefore, the "trumpet" regime 
holds, provided that 



1/N V «£</<!, 



(24) 



which means that then one deals with weak adsorption and weak Dulling force. 



2.2 Dependence of the restoring force cp on the adsorption energy e s 

One can write the restoring force in scaling form, ip = - s — J-(e s ), with e s = Eg/ksT and JF{x) - a 
scaling function. Force ip stays constant in the course of the desorption process as long as at least 
one adsorption blob is still located on the attractive surface. The strength of ip can be determined 
from the plateau in the deformation curve: force vs. chain end position [23J . 

In an earlier work [21], we demonstrated that the fugacity per adsorbed monomer z*(e s ), which 
determines its chemical potential /i at j s = In z*, can be found from the basic equation 

§(a,H3Z*)Q(\,ii2Wz*) = l. (25) 

In Eq. (125]) a = <ft — I (<ft is the so-called crossover exponent which, according to the different 
simulation methods, ranges between ~ 0.48 and <p ~ 0.59; see the corresponding discussion in 
ref. [24J ) A = 1 — 7,2=2 — —0.343 (where 7^=2 is a universal exponent which governs the two- 
dimensional polymer statistics), w = exp(e s ), /12, and /X3 are so-called connective constants in 
two- and three-dimensional spaces respectively, and the polylog function $ is defined by the series 
$(cv, y) = Yurv=i y n / na - The values of /i 2 , /x 3 and 7^=2 could be found in the text-book [25] . 

In local equilibrium /i a d s should be equal to the chemical potential of detached monomer /idet = 
1 dF stv 



, where F stT is the free energy of the driven polymer chain so that /i a d s = /idet- The 
kbT oN 

expression for F str depends on the polymer model. The moving blob domain, shown in Fig. [TJ 
exerts on the last adsorbed monomer a force which is equal in magnitude and opposite in sign to 
the restoring force ip. Therefore, the free energy of a stretched polymer within the bead-spring (BS) 
model is given by [24] 

/ \ 1/" 
F str = -N(-^\ -JV1b/x 3 (26) 

On the other hand, within the freely jointed bond vector (FJBV) model, the partition function, 
corresponding to F str , reads [26] 



7 ( A sN ( k B T 

Z N = (47T 
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(27) 



In Eq. ( |27p the free chain partition function Z^ = (4tt) for ip — Y 0. A correspondence with the 
BS-model may be established bv the substitution 4ii — > u s in En. (1271) so that the FJBV- free enerev 



becomes 
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Taking into account Eqs. (1231) and fl28|) . one obtains eventually 

^ _ / [~ MA^Ce,))]" , BS model 

T [ft-H-lnt/isZ*^))), FJBV model l J 

where 7tL _1 (x) is the inversed of the function 1Z(x) = In [/i 3 sinh(a;)/x] and fi 3 z*(e s ) < 1. 

Below we consider separately the restoring force, Eq. (1291) . in the limits of strong and weak 

adsorption. 

2.2.1 Strong adsorption 

In the case of the strong adsorption (e a ^> 1), the solution of Eq. (1231) reads z*(e s ) — Ai^ 1 e~ <Es [24J, 
and for the BS-model one gets aip/ksT ~ [e s + ln^/zis)]". 

For the FJBV model at e s ^> 1 one has 1Z{x) ~ In (/x 3 e x /x) ~ x + ln/i 3 . From Eq. (1291) it then 
follows aip/ksT + In /x 3 = e s + ln/i 2 , or aip/ksT ~ e s + ln(/i 2 //i 3 . The final result for the restoring 
force in the strong adsorption limit takes on the form 

~ _ j [ £s + ^(^a/jUs)]" , BS model 
^ [e s + ln(/x 2 //x 3 ), FJBV model 

It is well known that under strong forces (and correspondingly strong adsorption) the FJBV-model 
is better suited for description of the simulation experiment results (see, e.g., Sec. 5 in our paper 
). That is why in Eq. [30] the second line is appears as the more appropriate result. 



2.2.2 Weak adsorption 

Close to the critical adsorption point, the solution of Eq. ( 1231) takes on the form 

z*(e s ) ~ /Uj 1 [1 - Ci(e, - t c ) l/4 '] [21], where C x is a constant. Then, Eq. (J2SD gives (up/k B T ~ 

Ci(e, - e c )^ for the BS-model. 

In the case of weak adsorption aip/ksT is small and 7^(x) ~ In [7x3(1 + x 2 /6)} ~ ln/x 3 + x 2 /6. 

Eq. (1291) then yields atp/ksT ~ -Bi(e s — e c ) 1//2(?i , where i?i is a constant. The unified result in the 

weak adsorption limit is 

~ = fc 1 {e a -e c ) v/ +, BS model 

^ I Si fe. - e.) 1 /^. FJBV model l j 



In fact, Eq. (1311) may also be derived differently. Namely, as already mentioned, the size of an 
adsorption blob is given by D = aC\{t s — t c )~ u ^ whereas the size of a tensile blob £(x = D) = 
k B T/(p. The condition for detachment may be set as £(x = D) = D which leads again to Eq. 
( l3"Tj) (whereby the dimensionless notation tp = cnp/k B T has been used). On the other hand, the 
relationship between D and tp, written as D(p/k B T ~ 1, leads to the adsorption blob size for the 
FJBV- model: D ~ (e s — e c ) _1 ' 2?i (where the second line in Eq. (13~T1) has been used). 

2.3 Stem and stem-trumpet regimes 

For strong, 1 < a(p/k B T < af/k B T, and moderate, aip/k B T < 1 < af/k B T, detachment force / 
one has to consider the "stem" (see Fig. [Tb) and "stem-trumpet" (see Fig. [lb) regimes, respectively. 

2.3.1 Stem regime 

In this case the moving domain is build up from a "stem" of the length M(t) (total number of 
desorbed monomers). The "stem" velocity is dM(t)/dt (with an accuracy of the geometrical factor 
C = cosO which is not important for scaling predictions), so that the balance of driving and drag 
forces yields 

, m <m =f . v (32) 

The solution of Eq. (|32j) in the dimensionless units reads 



M(t) = yj2(f-lp)t (33) 



Evidently, this takes on the same form as Eq. ( )23|) . M(t) = N\/t/(r d ), whereby the mean detach- 
ment (dimensionless) time (r d ) = N 2 /[2(f — £>)]. 

2.3.2 Stem-trumpet regime 

The "stem-trumpet" shape of the driven polymer is shown in Fig. [Tb. In this case the local density 
in the "trumpet" part p(x, t) is governed by the same Eq. (j7]) but with a boundary condition in the 
junction point (i.e., in the point where the "stem" goes over into the "trumpet" part). There, the 
blob size is £(x = R(t) — S(t),t) = a. Thus, one gets 

d(x = R(t)-S(t)A) = l (34) 



The boundary condition at x = D has the same form, Eq. (151) , as before. The solution again has 
the scaling form p(x,t) = W(u). The moving end now follows the law R(t) = D + S(t) + Udvt, 
where the prefactor u d can be fixed by using the condition W{u = u d ) = 1. Consequently, 



u d = v/^~ 2 )/ 2 [1 - ^ 1 "")/"] 



(35) 



In order to calculate the length of the "stem" portion, S(t), one can use the results of the 
previous subsection. The only difference lies in the fact that now the "stem - trumpet" junction 
point moves itself with velocity Ud/(2vt) and the driving force is / — 1. Therefore, the balance of 
driving and drag forces leads to the following equation for S(t) 



2Vt dt 



(36) 



The solution of Eq. fl36l) reads 

S(t) 



§ + 2(7-1)-^ 

4 v 2 



t 



(37) 



The resulting equation for the average height of the end monomer becomes 

u d 



R(t) = D + 



u 



-£ + 2(f-l) + — 

4 v ' 2 



t 



(3? 



The number of desorbed monomer can be easily calculated as an integral over the the number 
of monomers in the "trumpet" part and in the "stem" portion, i.e., 

M(t) = 



t / W{u)du + S{t) 
o 

Ud 



U 



f +2(/-l)-yl I 



{l-v)/v 



<r 



(39) 



where ua is given by Eq. ( 1351) . Evidently, all characteristic scales as well as the number of desorbed 
monomers vary in time as \/t. This shows that the average desorption time (r) is proportional to 
iV 2 as before. 

2.4 Dynamic Regimes 

The foregoing theoretical analysis was based essentially on the notion of tensile (Pincus) blobs as 
well as on the "auasi-static" approximation, i.e., on the eauating of driving and drag forces locallv 



(in space and time). Depending on the pulling / and restoring p> force strengths, one distinguishes 
a "trumpet" {l/N v < p> < J < 1), "stem-trumpet" (<p < 1 < / ), and a "stem" (1 < <p < J) 
regimes (see Fig. [2]). Notably, in all regimes the average number of desorbed monomers M(t) and 
the height of the first monomer R(t) evolve as a yt-txme universal law. Moreover, the evolution 
equation which governs the polymer beads density p(x,t) in the tensile blobs is analogous to the 
nonlinear porous medium equation which is commonly used for investigations of gas and fluid 
infiltration in the geological porous strata [TTJ [18]. Following this line we have demonstrated that 
the scaling analysis of such equation, based on symmetry of self-similarity, immediately brings about 
a v^-time universal law. 




1/N 



Figure 2: Different dynamic regimes of forced polymer chain desorption in terms of pulling / and 
restoring p> forces. The relation between <p and the adsorption energy e s is given by Eq. (1301) and 
Eq. (I3T1) for the strong and weak adsorption, respectively. The white strip along the diagonal 
indicates that in all cases the condition / > p is important. 



The value of p> depends on the adsorption energy e s and is described by different expressions, 
depending on the particular polymer chain model: bead-spring (BC), or freely jointed bond vectors 
(FJBV) models. For large adsorption energy (e s ^> 1), the restoring force is also large, leading to 
the "stem" regime. Close to the the critical adsorption point (for vanishing adhesion), the restoring 
force could attain a very small value, leading to a 'trumpet", or "stem-trumpet" regimes. In this 
case the role of fluctuations, which was neglected in our theoretical investigation, becomes important 
and could result in correction to the basic -y/t-time law. 

In the next section we will give results of our MC- and MD-simulation studv that were carried 



out in some of the foregoing forced-desorption (mainly "stem") regimes. 



3 Computational Models 
3.1 Monte-Carlo 



We have used a coarse grained off-lattice bead spring model to describe the polymer chains. Our 
system consists of a single chain tethered at one end to a flat structureless surface while the external 
pulling force is applied to other chain end. In the MC simulation the surface attraction of the 
monomers is described by a square well potential U w (x) = —e s for x < 5 = 0.125 and U w (x) = 
otherwise. Here the strength e s /kBT is varied from 0.6 to 3.6. The effective bonded interaction is 
described by the FENE (finitely extensible nonlinear elastic) potential. 
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(40) 



The non-bonded interactions between the monomers are described by the Morse potential. 
U M (r) 



£m 



exp(-2a(r - r min )) - 2exp(-a(r - r min )) 



(41) 



with a = 24, r min 



0.8, tM/ksT = 1. This model is currently well established and has been used 
due to its reliability and efficiency in a number of studies [27] of polymer behavior on surfaces. 





Figure 3: Snapshots from a MC simulation of chain with N = 128, adsorbed on a solid plane with 
energy of adsorption e s /k B T = 3.0. (a) at time t — before, and (b) - after a pulling force / = 3.0 
is applied perpendicular to substrate. 



We employ periodic boundary conditions in the y — z directions and impenetrable walls in the z 
direction. The lengths of the studied polymer chains are typically 64, 128, and 256. The size of the 
simulation box was chosen appropriately to the chain length, so for example, for a chain length of 
128, the box size was 256 x 256 x 256 . All simulations were carried out for constant force. A force 
/ was applied to the last monomer in the ^-direction, i.e., perpendicular to the adsorbing surface 

The standard Metropolis algorithm was employed to govern the moves with self avoidance au- 
tomatically incorporated in the potentials. In each Monte Carlo update, a monomer was chosen at 
random and a random displacement was attempted with Ax, Ay, Az chosen uniformly from the in- 
terval —0.5 < Ax, Ay, Az < 0.5. The transition probability for the attempted move was calculated 
from the change AU of the potential energies before and after the move as W = exp(— AtT/A^T). 
As for standard Metropolis algorithm, the attempted move was accepted, if W exceeds a random 
number uniformly distributed in the interval [0, 1]. As a rule, the polymer chains have been orig- 
inally equilibrated in the MC method for a period of about 10 6 MCS (depending on degree of 
adsorption e s and chain length N this period is varied) whereupon one performs 200 measurement 
runs, each of length 8 x 10 6 MCS. Various properties of the chain are then sampled during the 
course of the run. The starting configuration is replaced by a new sequence in the beginning of the 
next run. Two typical snapshots of chain configurations, immediately before the detachment starts, 
and some time later, are presented in Fig. |3j 

3.2 Molecular Dynamics 

In our MD-simulations we use a coarse-grained model of a polymer chain of N beads connected by 
finitely extendable elastic bonds. The bonded interactions in the chain is described by the frequently 
used Kremer-Grest potential, V KG (r) = V FENE (r) + V WCA (r), with FENE potential, Eq. Q3QJ, and 
a non-bonded repulsion term taken as Weeks- Chandler- Anderson (WCA) (i.e., the shifted and 
truncated repulsive branch of the Lennard- Jones potential) given by 

V WCA (r) = 4e [(a/r) 12 - (a/rf + 1/4] 6(2 1/6 a - r) (42) 

with 9{x) = or 1 for x < or x > 0, and e = 1, a = 1. The potential V KG (r) has got a minimum 
at bond length r bond us 0.96. 



The substrate in the present investigation is considered simply as a structureless adsorbing 
plane, with a Lennard- Jones potential acting with strength e s in the perpendicular x-direction, 
V w {x) = 4e s [(a/x) 12 - [a/xf]. 

The dynamics of the chain is obtain by solving a Langevin equation of motion for the position 
r n = [%n, Vni z n ] of each bead in the chain, 

mt n = ¥{ + F^ CA - 7*« + Rn(t) + Wntf, (jl, . . . , N) (43) 

which describes the Brownian motion of a set of bonded particles whereby the last of them is 
subjected to external (constant) stretching force f = [/, 0,0]. The influence of solvent is split into 
slowly evolving viscous force and rapidly fluctuating stochastic force. The random, Gaussian force 
R n is related to friction coefficient 7 by the fluctuation-dissipation theorem. The integration step is 



0.002 time units (t.u.) and time in measured in units of ymcr 2 /e, where m denotes the mass of the 
beads, m — 1. The ratio of the inertial forces over the friction forces in Eq. (1431) is characterized 
by the Reynolds number Re = ^Jmej^a which in our simulation falls in the range Re = 0.1 -j- 4. In 
the course of simulation velocity- Verlet algorithm is used to integrate equations of motion ( 1431 . 

4 Simulation results 

In order to check the validity of our theoretical predictions, we carried out extensive computer 
simulations by means of both Monte-Carlo (MC) and Molecular Dynamics (MD) so that a com- 
parison can be made between the overdamped dynamics of a polymer (MC) and the behavior of 
an inertial (underdamped) system (MD). For ultimate consistency of the data, we performed also 
MD simulations with very large friction coefficient 7 = 10 whereby the obtained data was found to 
match that from the MC computer experiment. 

First of all, it would be interesting to verify that trumpet-like structures, which have been dis- 
cussed thoroughly in SecJ21 can be seen in our MD-simulation experiment. Figured] represents such a 
visualization where the beads density distributions are given (as a function of ^-coordinate as well as 
the radial coordinate calculated with respect to the first monomer, i.e. r = 1/(2/1 — y) 2 + (z 1 — z) 2 ) 
for the weak , e s /ksT = 1.9, (left panel) and strong, e s /ksT = 5, (right panel) adsorptions. Each 
densitv plot was implemented as a result of averages over many runs under the fixed fraction of 
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Figure 4: Polymer density distributions of beads plotted for different adsorption regimes: (a)-(c) 
E s /{kBT) = 1.9, (d)-(f) e s /{k,BT) = 5.0 and various percentage of detached monomers (10%, 50% 
and 90%). The a;-axis coincides with direction of pulling force whereas r = y/(y\ — y n ) 2 + {z\ — z n ) 2 , 
(n — 1, . . . , N) is the radial component of beads calculated with respect to the first monomer. Here 
N = 100, T = 0.1 and 7 = 0.25. 



desorbed monomers: 10% (upper row), 50% (middle row) and 90% (bottom row). As one can 
see, the moving domain is developed into the trumpet-like configuration in both, weak and strong 
adsorption, cases. Moreover, the "trumpet" is tilted as it was argued in Sec. 12.1.11 due to its lateral 
velocity and the corresponding Stokes friction. An important point is that the tilt angle stays 
constant which substantially simplifies the theoretical consideration in Sec. 12.1.11 Recall that the 
density in the blob n&(x) and the blob size £(x) are related as n^x) oc l/[£,(x)] 3 ~ 1 / u oc l/[£(x)] L3 , 
i.e., the larger the density is, the smaller blob size. This tendency can be seen in Fig. H] where the 
most dense portion (red color) comes close to the tip of the "trumpet" as it should be. Moreover, 



for the strong adsorption case (see Fig. SU-f ) this most dense portion becomes more extended which 
might mean a "stem-trumpet" formation. As a result, the density plots, given in Fig. HJ provide a 
good evidence for the existence of the trumpet-like structures, suggested and discussed in Sec. El 
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Figure 5: (a) MD data for the mean desorption time (r^) of a chain vs chain length N, plotted 
for different strength of pulling force /, adhesion strength eg/ksT, and friction 7. Symbols denote 
the results of simulation while solid lines represent fitting curves (r^) oc N 13 which were found in 
the range N = 50 -r 250. The value of the exponent /3 depends on the dynamics: underdamped 
f3 ~ 1.75 ± 0.1 (full symbols) and overdamped 1.96 ± 0.03 (empty symbols), (b) Variation of (tA 
with iV at different strength of adhesion e s /ksT from the MC data. Note that here the critical 
adsorption energy e c /kBT m 2.0. In the strong adsorption regime (3 = 2.0, whereas close to the 
critical point (3 = 1.83. 



As one of our principal results, it was found that the chain length dependence of the mean 
detachment time (rA follows a scaling law (tA ~ N", Fig. [5j where /3, the detachment exponent, 
depends on the type of dynamics (underdamped or overdamped) as well as on the adsorption 
strength (i.e., on the proximity of e s to the adsorption critical point, e c ). Fig. [5h shows result from 
the MD-simulation for different values of e s , the pulling force /, and the friction coefficient 7. For the 
underdamped dynamics, one finds a detachment exponent (3 = 1.75=1=0.1 whereas for the overdamped 
one (3 = 1.96 ± 0.03. The last value is very close to the theoretical prediction which immediately 
follows from Eq. (133|) . (r^) oc T N 2 /(f — <p), i.e., (3 = 2. Recall that this theoretical prediction 
neglects any effects of inertia and describes overdamped dynamics. In addition, detachments are 
considered in the "stem" scenario, that is, for relatively large adsorption energy and pulling force. 

This conclusion is supported by the MC-simulation results (which represent the overdamped 



limit) shown in Fig. \5jp. For the relatively strong adsorption one finds (3 = 2.02 ± 0.01, which is 
very close to the theoretical prediction. On the other hand, in the vicinity of the adsorption critical 
point the detachment exponent appears to be smaller (3 = 1.83 ± 0.01. Indeed, one may suggest 
that the observed decrease is due to fluctuations which become stronger close to e c . As mentioned 
above, fluctuations are not taken into account in our theoretical treatment which is based on the 
"quasistatic" approximation. As a matter of fact, in the cases of "trumpet" and "stem-trumpet" 
scenarios (relatively weak adsorption energy and pulling force), the fluctuations alter the value of 
the detachment exponent so that /3 < 2. It is, however, conceivable too that (3 = 1.83 may also be 
due to finite-size effects in the vicinity of the critical adsorption point where the size of adsorption 
blobs becomes large. 
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Figure 6: MD data on the inverse mean desorption time 1/(t^) of a chain with N = 50 plotted as a 
function of: a) pulling force /, and b) adsorption strength e s /kBT. Here 7 = 0.25. c) Inverse mean 
detachment time N 2 / {t^) (scaled with N 2 ) against e s from MC data with / = 3.0 and e s /ksT = 3.0. 



In order to check the /— and e s — dependence of (r d ), which the relation (r d ) oc r N 2 /(f — <p) 
predicts for the the "stem" regime, we have plotted l/fa) vs. /, Fig. EK, and vs. e s , Fig. [6b. In 
Fig. [6k one can see a perfect straight line as expected. As for the e s -dependence, we recall that 
according to Eq. (!30j) . in the strong adsorption limit the restoring force (p ~ c\(e a + c 2 ) where c\ and 
C2 are some constants. As can be seen from Fig. [6b, in the MD-simulation this holds for e s > 5. On 
the other hand, in the MC simulation, the strong adsorption regime starts at e s > 2.5 as is evident 
from Fig. [6b. 

Fig. [7] demonstrates the validity of the y/i-law for the time evolution of the mean number 
of desorbed monomers, M(t), given by Eq. (1331) . According to this relationship, the quantity 
jM 2 (t)/(f — tp) should be a linear function of time. 

In order to fix the value of the restoring force id, we have carried out an independent simulation 
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Figure 7: (a) Mean squared number of desorbed monomers 7M 2 (£)/[(/ — if)} for a chain composed 
of N = 100 particles shown as a function of time for different sets of / and 7 - MD data. Here 
Eg/hsT = 5.0. The inset displays the (equilibrium) relationship between pulling force / and the 
resulting mean height (H) of the first monomer in a chain with N = 50. The estimated value of the 
plateau height (which is the restoring force in the course of detachment) is ip = 0.44. (b) Variation 
of M(t) with (dimensionless) elapsed time (measured in units of t/(Td)) for different polymer lengths 
N from MC data. Here e s /k B T = 3.0. Dashed line denotes the expected \7i-behavior, cf. Eq. ( 1331) . 



to find how the mean height of the first monomer (H) depends in equilibrium on the pulling force 
/ for a previously adsorbed (end-grafted) chain. The result of this is shown as an inset in Fig. 
[7Ja where the value of plateau height yields the restoring force ip = 0.44 at adsorption energy 
e s = 5.0. By making use of this value we were able to superimpose three curves for different forces, 
/ = 2.5, 5, 10 in the overdamped case (7 = 10.0) on a single master curve as shown in Fig. [7_k 
Evidently, in the case of underdamped dynamics this collapse on a single master curve does not 
work. This deviation from the theoretical predictions indicates again that they apply perfectly 
only to the case of overdamped dynamics. Finally, the respective M(t) vs. time relationship from 
the MC-simulation is shown in Fig. [7b. One should thereby note that in the starting equilibrium 
configuration of the adsorbed chain about 24% of the monomers reside in loops so they are actually 
desorbed even before the pulling force is applied. Evidently, as in the MD data for very large 
friction, the curves for different chain lengths collapse on a single master \/t-curve in agreement 
with Eq. flM}. 



5 Plausible experimental validation 

Now we are in a position to discuss how the detachment kinetics considered in this paper could 
be validated in a laboratory experiment. To this end we recall that using conventional AFM, 
optic (OT), or magnetic tweezers (MT), the tracking is performed by practically macroscopic, of 
micrometers large, agents: by cantilever in AFM, or by e.g. polystyrene beads in OT and MT 
cases. On the other hand, the trumpet-like structures are the result of delicate interplay between 
pulling and drag (or friction) forces. Obviously, the Stokes friction of the cantilever as that of the 
polystyrene beads exceeds largely the friction of the tested polymer chain itself and thus suppresses 
the observation of the predicted time-behavior. One could solve this problem only if the the tracking 
agent is of the size of a single monomer, e.g., when it is a charged ionic end-group like a carboxylate 
group placed in an external electric field. Moreover, the polymer end should be fluorescently labeled 
in order to measure the polymer end spatial location. 

One way to do this is by means of total internal reflection fluorescence microscopy (TIRFM) 
which is a powerful technique aimed at imaging fluorescent species close to an interface [28J. A 
possible schematic setup is shown in Fig. E] 
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Figure 8: Schematic representation of micromanipulation (detachment) experiment using TIRF 
microscope. 



First of all, let us estimate the required tracking forces. In the strong adsorption limit aip/ksT 



0.5. The Boltzmann constant ks = 1.4 • 10~ 23 J/K, the temperature T = 300 K, a bond length 
a = 10~ 10 m so as a result p = 2.1-10 -11 N = 21 pN. In the case of weak adsorption, ip ~ (ksT /aN v ), 
that is, for N = 10 3 one has ip ~ 8.4 ■ 10~ 13 N « 1 pN. As a result, one deals typically with forces 
1 pN < p < 20 pN. 

Now, one could estimate which electric field E should be applied to electrodes, of the capacitor 
shown in Fig. so as to produce a tracing force 1 pN. The tracking force is / = qE, where q is the 
charge of the ionic end-group. The typical ion charge q = 1.6 • 10~ 19 C , so that for E = 5 • 10 6 V/m 
one gets / = 5-1.6-10~ 19 -10 6 ~ 10 -12 N = 1 pN. This value of the electrical field is much smaller than 
the dielectric strength of benzene, E st = 165 • 10 6 V/m, or even of distilled water E st = 70 ■ 10 6 V/m 
[29] . i.e., these fluids could be used as dielectrics between the plates of the capacitor. We underline 
that a totally stretched polymer of chain length, say N = 10 3 , has a characteristic size of hundred 
nanometers, so that the distance between plates of the capacitor should not be larger than l-=-10/im. 
This means that in order to create the electric field E = 5-10 6 V/m one should apply to the capacitor 
a voltage 5 -j- 50 V, i.e., no high- voltage facility are necessary. 

From a chemical point of view, the sample represents a substrate with end-grafted polymers 
at a low graft density (i.e., the chains are in a "mushroom" conformation). The synthesis of, e.g., 
polystyrene end-grafted chains could be performed by surface-initiated radical polymerization [30J, 
followed by the labeling of the other chain end through carboxyl groups and fluorescent dyes. These 
estimations show that an electrostatic tracking experiment provides a way to test and verify our 
main theoretical results. 

6 Conclusion 

In this work we consider the force-induced desorption of a polymer from a structureless adhesive 
surface both analytically and by means of two distinct methods for computer simulation. 

We have shown theoretically that there are three dynamical regimes of a polymer chain detach- 
ment. Depending on the pulling / and restoring ip forces, one can discriminate between a "trumpet" 
(1/N V <C (p < f < 1), "stem-trumpet" (<p < 1 < /), and "stem" (1 < <p < f) regimes of desorption. 

Remarkably, in all these cases the time dependence of the number of desorbed monomers M (t) 



and the height of the first monomer (i.e., the monomer which experiences the applied external 
pulling force) R(t) follows an universal v^Taw (even though this is not a diffusion phenomenon). 
There is, however, a common physical background with the well-known Lucas- Washburn y^-law of 
capillary filling J2I] as with the ejection kinetics of a polymer chain from a cavity (virus capsid) 
|32j . In these seemingly different phenomena there is always a constant driving force (meniscus 
curvature, or polymer entropy) which acts against a gradually changing drag force (friction) in the 
course of the process. 

We discovered an interesting similarity between the differential equation governing the monomer 
density variation p(x, t) in time and space, and the nonlinear porous medium equation (PME) 
[T7] [18] . This makes it possible to use the self- similarity property of PME and derive rigorously the 
-y/t-law. 

Our extensive MD- and MC-simulations of the detachment kinetics support this finding as well 
as the {jd) oc N 13 scaling with (3 = 2 of mean detachment time r^ with chain length N. The 
theoretically predicted dependence of (r^) on pulling force / and adsorption energy e s appears in 
perfect agreement with the simulation results. As noted above, the consistency between theory and 
computer experiment is well manifested in the case of overdamped dynamics and strong adsorption. 
Moreover, by means of MD-simulation we have shown that the beads density distribution plots 
support the notion of trumpet-like tilted tensile blob structures. One can envisage not only tilted 
but, e.g., bended (horn-like) structures where the local lateral velocity depends on the x-coordinate. 
We leave this more complicated case for future investigation. 

The deviations in the exponent /3 due to inertial effects in the underdamped dynamical regime 
(/3 ~ 1.75 found in the MD simulation) as well as close to the adsorption critical point e c ((3 = 1.83, 
in the MC simulation) are also challenging. It appears possible that inertial effects and fluctuations 
close to the CAP (both neglected in our theoretical treatment) may lead to similar consequences 
for the dynamics of polymer desorption, manifested by the observed decrease in the detachment 
exponent /3. In both cases this suggests that polymer detachment is facilitated. This is clearly a 
motivation for further studies of the problem. 



Acknowledgements 

We thank K.L. Sebastian, H.-J. Butt, K. Koynov, and M. Baumgarten for helpful discussions. 
A. Milchev is indebted to the Max-Planck Institute for Polymer Research in Mainz, Germany, 
for hospitality during his visit and to CECAM-Mainz for financial support. This work has been 
supported by the Deutsche Forschungsgemeinschaft (DFG), grant No. SFB 625/B4. 

References 



[i] 

[2] 
[3] 
[4] 

[5] 
[6] 

[7 



[9 
[10 

[11 

[12 
[13 



Ritort, F. J. Phys. Cond. Mat. 2006, 18, R531. 

Bhattacharjee, S. M. J. Phys. A 2000, 33, L423. 

Marenduzzo, D.; Trovato, A.; Maritan, A. Phys. Rev. E 2001, 64, 031901. 

Orlandini, E.; Bhattacharjee, S.; Mernduzzo, D.; Maritan, A.; Seno, F. J. Phys. A 2001, 34, 
L751. 

Kierfeld, J. Phys. Rev. Lett. 2006, 97, 058302. 

Sebastian, K.L. Phys. Rev. E 2000, 62, 1128. 

Marenduzzo, D.; Bhattacharjee, S.M; Maritan, A.; Orlandini, E.; Seno, F. Phys. Rev. Lett. 
2002, 88, 028102. 

Bhattacharya, S.; Rostiashvili, V. G.; Milchev, A.; Vilgis T. Phys. Rev. £2009, 79, 030802(R). 

Bhattacharya, S.; Rostiashvili, V. G.; Milchev, A.; Vilgis, T. Macromolecules 2009, 4®, 2236. 

Bhattacharya, S.; Milchev, A.; Rostiashvili, V. G.; Vilgis, T. Eur. Phys. J. E 2009, 29, 285. 

Skvortsov A. M.; Klushin L. I.; Fleer G. J.; Leermakers, F. A. M. J. Chem. Phys. 2010, 132, 
064110. 

Brochard-Wyart, F. Europhys. Lett. 1993, 23, 105. 

Sakaue, T. Phys. Rev. E 2007, 76, 021803. 



[14] Sakaue, T. Phys. Rev. E 2010, 81, 041808. 

[15] Sebastian, K.L.; Rostiashvili, V.G.; Vilgis, T.A.; EPL 2011, 95, 48006. 

[16] Simsen, J.; Gentile, C.B. Nonlinear Anal. 2009, 71, 4609. 

[17] Vazquez, J.L. The Porous Medium Equation, Clarendon Press, Oxford, 2007. 

[18] Barenblatt, G.I.; Entov, V.M.; Ryzhik, V.M. Theory of Fluid Flows Through Natural Rocks, 
Kluwer Academic Publishers, London, 1990. 

[19] de Gennes, P.G.; Pincus, P. J. Phys, (France) Lett. 1983, 44, L241. 

[20] de Gennes, P.G. Scaling Concept in Polymer Physics, Cornell University Press, Ithaca, 1979. 

[21] Sakaue, T.; Yoshinaga, N. Phys. Rev. Lett. 2009, 102, 148302 

[22] Serr, A.; Netz, R. Europhys. Lett. 2006, 73, 292. 

[23] Bhattacharya, S.; Milchev, A.; Rostiashvili, V.G.; Vilgis, T.A. Eur. Phys. J. E 2009, 29, 285. 

[24] Bhattacharya, S.; Milchev, A.; Rostiashvili, V.G.; Vilgis, T.A. Macromolecules 2009, 42, 2236. 

[25] Vanderzande, C. Lattice Model of Polymers, Cambridge University Press, Cambridge, 1998. 

[26] Sheng, Y.-J.; Lai P.-Y. Phys. Rev. E 1997, 56, 1900. 

[27] Milchev, A.; Binder, K. Macromolecules 1996 29, 343. 

[28] Mashanov G.I.; Tacon D.; Knight A.E.; Peckham M.; Molloy J. E. Methods, 2003, 29, 142. 

[29] CRC Handbook of Chemistry and Physics, 93rd Edition, 2012-2013. 

[30] Y. Tsujii, K. Ohno, S. Yamamoto, A. Goto, T. Fukuda, Adv. Polym. Sci. 2006, 197, 1. 

[31] Dimitrov, D. L; Klushin, L.; Milchev, A.; Binder, K. Phys. Fluids2008, 20, 092102. 

[32] Milchev, A.; Klushin, L.; Skvortsov, A.; Binder K. Macromolecules 2010, 43, 6877. 



Polymer Detachment from Adsorbing Surface: Theory, 
Simulation and Similarity to Infiltration into Porous Medium 

by Jaroslaw Paturej, Andrey Milchev, Vakhtang G. Rostiashvili, and Thomas A. Vilgis 
TOC graph: 




Figure 9: A snapshot from MC simulation of a partially desorbed polymer chain with N = 128 
monomers under pulling force / = 3.0 and surface adhesion e^/fe^T = 2.50. 



